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Catalytically active colloids maintain non-equilibrium conditions in which they produce and de¬ 
plete chemicals and hence effectively act as sources and sinks of molecules. While individual colloids 
that are symmetrically coated do not exhibit any form of dynamical activity, the concentration 
fields resulting from their chemical activity decay as 1/r and produce gradients that attract or 
repel other colloids depending on their surface chemistry and ambient variables. This results in a 
non-equilibrium analogue of ionic systems, but with the remarkable novel feature of action-reaction 
symmetry breaking. We study solutions of such chemically active colloids in dilute conditions when 
they join up to form molecules via generalized ionic bonds, and discuss how we can achieve structures 
with time dependent functionality. In particular, we study a molecule that adopts a spontaneous 
oscillatory pattern of conformations, and another that exhibits a run-and-tumble dynamics similar 
to bacteria. Our study shows that catalytically active colloids could be used for designing self- 
assembled structures that posses dynamical functionalities that are determined by their prescribed 
3D structures, a strategy that follows the design principle of proteins. 


I. INTRODUCTION 

It is fascinating to learn how Nature exploits sophisti¬ 
cated mechanisms to keep time [T] , since having a system 
with temporal structure such as spontaneous oscillations 
in an over-damped inertia-less world is far from trivial. 
Oscillators are ubiquitous in biology [2], and they often 
play a vital role, such as the example of mitotic spindle 
oscillations that regulate the key function of chromosome 
separation during cell division [3]. Oscillatory behaviour 
typically arises through coordination of many stochas¬ 
tic components [4], and could have remarkable charac¬ 
teristics, such as self-tuning to criticality in the case of 
hair-bundles m- Another notable example of temporal 
structure is the run-and-tumble behaviour in the swim¬ 
ming pattern of microorganisms: the trajectories have 
relatively long run segments that are intercalated with 
burstlike tumble events when the orientation of the swim¬ 
ming microorganism is completely randomized [8]. This 
behaviour has been extensively studied for many bacteria 
that are propelled by rotary motors and stiff flagella, but 
has also been observed for eukaryotic alga that swim via 
synchronized beating of flexible flagella [9l [To] . For E. 
coli, the coupling between a highly sensitive chemotactic 
circuitry and the motility mechanism has been unraveled 
to exquisite detail m- 

It will be desirable to make artificial microscopic de¬ 
vices with autonomous dynamic functionality, such as 
those described above, through self-assembly of nano¬ 
scale building blocks, as is the case with the biologi¬ 
cal examples. Since the turn of the century, there have 
been numerous manifestations of synthetic devices with 
mechanical functionality using a number of different ap¬ 
proaches, including self-assembled DNA nano-structures 
mini, actuators that are triggered by global oscillatory 


chemical reactions [HUS], and magnetically actuated ar¬ 
tificial micro-swimmers Hi nil and artificial cilia [T8ti2Q] . 
A particular class of such active systems [21] that uses 
long-range phoretic interactions [22l [23] has been shown 
to lead to the emergence of collective activities including 
pattern formation and spontaneous oscillations [24l [25] , 
as well as swarming [26] , in the case of homogenous so¬ 
lutions. In heterogenous mixtures, chemotactic interac¬ 
tions could lead to spontaneous formation (self-assembly) 
of active colloidal molecules with nonequilibrium activity 
or function that is determined by their stable 3D struc¬ 
ture; a notion that follows the design principle of proteins 

m- 


The model for self-assembled active molecules pre¬ 
sented in Ref. m was used to demonstrate the forma¬ 
tion of low weight colloidal molecules that exhibit dif¬ 
ferent types of static activity, namely translational and 
rotational self-propulsion. The type of activity depends 
on the global parameters and the number of colloids that 
composed each assembly. Here we show that the model 
allows for the self-assembly of colloidal molecules that 
adopt time-dependent configurations, increasing the pos¬ 
sible states that active matter can manifest. Specifically, 
we show that a molecule can switch spontaneously be¬ 
tween active and passive states, mimicking a behaviour 
that resembles that of bacteria in their run-and-tumble 
motion. Also, we show that colloidal molecules can sus¬ 
tain spontaneous oscillations with frequencies that can 
be tuned by the global parameters, and hence act as mi¬ 
croscopic clocks. 
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II. THE MODEL 


tide 2 due to the activity of particle 1 is equal to 


We consider a suspension of spherical colloidal parti¬ 
cles that have a catalyst coating on their surfaces, in a 
solution of reactants that are catalytically converted into 
products at the surfaces of the colloids. We choose a mix¬ 
ture of particles with uniform coating, and, for simplicity, 
use a model in which the catalytic activities of the col¬ 
loids are simplified into net production or consumption 
of chemicals with given rates, which we denote as surface 
activity [28143^ . An isolated colloid with surface activ¬ 
ity a, which is positive (negative) when the net activity 
amounts to production (consumption) of chemicals, will 
produce a stationary spherically symmetric concentra¬ 
tion profile C around it given by C = Co + acr^/(4Dr), 
where a is the diameter of the colloid, D is the (nominal) 
diffusion coefficient of the chemicals, r is the distance to 
the centre of the colloid, and Co is a reference concen¬ 
tration at infinity. Since the Brownian diffusivity of the 
colloid, which we denote as Dc, is much smaller than D, 
the solute concentration profile relaxes very quickly to 
a comoving cloud when a colloidal particle moves. For 
simplicity, we ignore the effect of advection, which at fi¬ 
nite Peclet numbers will distort this cloud, because it 
will affect the self-propulsion velocity only within a pre¬ 
factor of order unity [33]. We also ignore the possibility 
of spontaneous symmetry breaking at large Peclet num¬ 
bers [34] and the anomalous super-diffusion at relatively 
short time scales [35] . 

In presence of a concentration gradient a coated 
particle will move with a velocity equal to V = 
— f dSVjjC, where /i is the phoretic mobility, which 
itself depends on the molecular interactions between the 
coated surface and the dissolved chemical. As in the case 
of a, this coefficient too can be both positive and nega¬ 
tive. Due to these phoretic effects, if two or more colloidal 
particles are placed close to each other, they will acquire 
drift velocities that can be regarded as effective nonequi- 
librium interactions. Remarkably, these interaction are 
asymmetric as the gradient each particle creates is con¬ 
trolled by its a whereas its response the gradient created 
by others is controlled by fi] two independent material 
parameters. More explicitly, the drift velocity of particle 
2 due to the activity of particle 1 will be proportional 
to 0 ^ 1 /i 2 whereas the drift velocity of particle 1 due to 
the activity of particle 2 will be proportional to 0 ^ 2 Mi • 
When the effective interactions between the particles are 
not symmetric, the system cannot reach an equilibrium 
state because the condition of detailed balance will not be 
fulfilled. This can manifest itself in the form of frustra¬ 
tion that leads to nonequilibrium fluxes. For example, 
if 0^1 >0 and /ii > 0? while 0^2 > 0 and M 2 < 0, 1 
will be repelled by 2 but 2 will be attracted by 1. In a 
mixture, we generically have (aiM 2 7 ^ <^iM 2 and, hence, 
nonequilibrium colloidal activity. For dissimilar colloids, 
the condition of detailed balanced can only be satisfied 
by fine-tuning. 

In the far-held approximation the drift velocity of par- 



aiiX2<y'^ ri2 
247rZ) |ri2|3’ ^ ^ 


where ri 2 = r 2 — ri. The effective interaction is similar 
to unscreened electrostatic interaction between charged 
colloidal particles in a huid. Hence, we can regard the 
two parameters as generalized charges: a is responsible 
for the production of the held and m controls the response 
to the held. When we are not in the far-held limit, the 
concentration prohle should be calculated by solving the 
dihusion equation with the appropriate boundary condi¬ 
tions, and the resulting drift velocities will be modihed 
in the near-held region and in general will not be pair¬ 
wise additive. However, these many-body and proxim¬ 
ity ehects will only quantitatively change the predictions 
that emerge from an analysis based on the calculation 
of the interactions in the far-held limit, where it is pair¬ 
wise additive, and not ahect the qualitative features. In 
particular, those ehects do not alter the main property 
of the colloidal interactions, namely, the asymmetry of 
action versus reaction. The far-held approximation of 
the Coulomb interactions augmented with short range 
steric repulsion between the particles make the basis of 
the highly successful and widely used restricted primitive 
model (RPM) for charged colloids [36]. Here we adapt 
this model to our nonequilibrium generalized Coulomb¬ 
like system. We also ignore hydrodynamic interactions 
between colloids as they respect the action-reaction sym¬ 
metry and appear as sub-dominant terms in the far-held 
limit. In the near held, the hydrodynamic interactions 
will modify the attraction or repulsion speeds between 
colloids due to lubrication ehects, changing only quanti¬ 
tatively the predictions made using this simplihed model. 

We thus perform a Brownian dynamics simulation of 
active colloidal particles by solving the following stochas¬ 
tic equations 


^=VoJ2 + E (2) 

where ^ is a random velocity represented by white noise 
of intensity 2Dc with Dc being the (passive) dihusion 
coefhcient of the colloidal particles, and Uik = —'^ki 
is the steric repulsion term that prevents particles from 
overlapping. Using ao and mo as characteristic values 
for the surface parameters, we have made the surface ac¬ 
tivity and mobility parameters dimensionless by dehning 
a = ajaQ and Ji = fijfiQ. This gives us an overall velocity 
scale of Vb = and the dimensionless noise intensity 

of Dc = DcjVoa^ which we regard also as a dimensionless 
temperature. For simplicity the interactions are com¬ 
puted assuming that the solute molecules could dihuse 
in three dimensions (with the ehective potential decay¬ 
ing as 1/r) but they colloids are restricted to move in 
two dimensions. To simulate the excluded volume term, 
we follow the strategy devised in Ref. m- We advance 
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the system in time steps of 6t = O.OOla/Vo, and at each 
time step, compute the drift velocity for each particle 
by summing all pairwise interactions. The particles are 
then advanced using a forward Euler scheme that adds 
the corresponding noise term. At this stage, we identify 
pairs of overlapping particles and reflect them, in order, 
by the same distance that they overlap. We repeat the 
procedure until there are no remaining overlaps. 

In what follows we consider a binary mixture, with 
both species—labelled A and B —having the same diam¬ 
eters, but different charges. We explore the case in which 
A and B particles are mutually attracted (although with 
different intensities), while identical particles repel each 
other. The characteristic values ao and fio are chosen 
such that for B particles Ab = Jib = —^ and /I^ > 0. 
Without loss of generality, we choose Sa > 1 (the oppo¬ 
site case is obtained by exchanging the role of A and B 
particles). 


III. OSCILLATORY INSTABILITY 

The absence of action-reaction symmetry in the effec¬ 
tive interaction between colloids implies that the system 
does not necessarily evolve toward minima of effective 
nonequilibrium free energies: it is in principle possible 
to have residual dynamics in the long time limit, e.g. in 
the form of limit cycles. The limit cycles will be asymp¬ 
totic time-dependent stable solutions where the system 
may be trapped (despite the Brownian noise) until large 
perturbations occur, possibly in the form of collisions of 
active molecule with other molecules in the present case. 

We are interested in finding low weight molecules that 
can sustain stable oscillatory motion. To isolate the oscil¬ 
lation from other kinds of activity, we consider molecules 
with polar and plane-reflection symmetries to avoid oc¬ 
currence of translational and rotational self-propulsion. 
Figure [^presents a sample configuration of the molecule 



FIG. 1. Configuration of a molecule with polar and plane- 
reflection symmetries, where the colloids are free to move with 
respect to one another subject to the constraint that A and 
B colloids remain in contact. It has Na = 4 and Nb = S 
colloids and there are 10 degrees of freedom. The oscillations 
of this molecule are characterized by the parameter y = ([^^i ~ 
1 * 2 ] X [vcm ~ 1 * 0 ]) • z, where Vcm is the position of the centre of 
mass, ri /2 are the positions of the two extreme particles, and 
1*0 = ( 1*1 + r 2 )/ 2 , as labelled in the figure. 


with the aforementioned symmetry, which we will con¬ 
sider. As A and B colloids mutually attract, they will 
remain in contact if the noise intensity is small. There¬ 
fore, although simulations are performed with full dy¬ 
namics, we can introduce a further simplification for the 
purpose of performing the analysis, and consider the case 
where the particles are restricted to remain in contact 
and use d’Alembert’s principle to derive the dynamics of 
the remaining coordinates (see Appendix 0 - To param¬ 
eterize the motion of this molecule we chose the following 
generalized coordinates: the x and y coordinates of the 
leftmost colloid and all the angles that give the position 
of the colloids relative to the previous one, going from left 
to right, always preserving the condition that A and B 
colloids are in contact. The equations of motion that re¬ 
sult after using d’Alembert’s principle are first solved to 
find the equilibrium configuration (regardless of it stabil¬ 
ity). Then, the equations are linearized around equilib¬ 
rium, which are analyzed to find relaxation modes with 
exponential time dependence The dimension¬ 

less eigenvalues A depend on the values of the charges, 
and we find that a Hopf bifurcation can take place. Fig¬ 
ure [^present the eigenvalues at the bifurcation, where a 
pair of complex eigenvalues cross the imaginary axis, and 
acquire positive real parts. There are three null eigenval¬ 
ues associated to the translational and rotational sym¬ 
metries. At the bifurcation all other eigenvalues are neg¬ 
ative, corresponding to damped motion. Consequently, 
as it is generically the case, in the neighborhood of the 
Hopf bifurcation the stable fixed point associated with 
the straight molecule gives birth to a small stable limit 
cycle. The limit cycle corresponds to a bending oscilla¬ 
tion of the molecule, which, in turn, generates a periodic 
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FIG. 2. Real and imaginary parts of the dimensionless eigen¬ 
values A of the linearized motion of the A^Bq molecule (shown 
in Fig. for a A = 1.5 and JIa — 0.7, where two complex 
conjugate eigenvalues cross the imaginary axis, with imagi¬ 
nary parts of ±0.127. Only the region close to the imaginary 
axis is shown. Not shown are the three remaining eigenval¬ 
ues, which are real and negative. The null eigenvalue has 
degeneracy three. Increasing a a or decreasing JIa moves the 
two critical eigenvalues to the positive real region, inducing 
an instability that is saturated with non-linear terms leading 
to the development of a limit cycle at the Hopf bifurcation. 
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FIG. 3. Configurations of the molecule that undergoes oscil¬ 
latory motion at different phases of the cycle, with a a = 1.5 
and JIa = 0.6. The full motion is presented in the supplemen¬ 
tary material [39] . 


oscillation of the center of mass in the transverse direc¬ 
tion. No rotation is obtained and the net translation in 
one cycle vanishes. Figure displays the configurations 
at different phases of the cycle and the full motion is pre¬ 
sented in the supplementary material [39] . The period of 
the oscillation is given by the imaginary part of the two 
critical eigenvalues and the amplitude of the oscillation 
increases with the distance to the critical point. A movie 
of the self-assembly process that leads to the formation of 
the molecule from an initial random dispersion of 

colloids is presented in the supplementary material [39] . 

The quality of the microscopic oscillator can be studied 
by analyzing the temporal evolution of a simple observ- 
able. We define x = ([ri - r 2 ] x [rcm - (ri + i’2)/2]) • z, 
where Vcm is the position of the centre of mass, and ri /2 
are the positions of the two extreme particles (see Fig. 
[^. This observable measures the instantaneous defor¬ 
mation of the cluster from the linear configuration, while 
being rotationally and translationally invariant. Figure 
[^top presents the typical time evolution of x and its 
power spectrum for different noise intensities. There is 
a clear peak at cj ^ 0.13Vb/cr, consistent with the fre¬ 
quency predicted by the eigenvalue analysis (Fig. [^. 
There is a second peak, which corresponds to the third 
harmonic (there is no second harmonic by symmetry). 
As expected, the width of the central peak decreases by 
decreasing the noise intensity. At Dc = b x 10 ^ the fre¬ 
quency is cj = (0.1326 ± 0.0002) Vb/cr, whereas for noise 
intensities larger than Dc = 2.5 x 10“^, the fluctuations 
become large and break the cluster, which no longer os¬ 
cillates. The oscillations also become sharper by moving 
away from the transition, as it is shown in Fig. [^bottom, 
where the power spectrum peak increases by increasing 

The colloidal oscillator ^ 4^3 is the smallest molecule 
we found where the first eigenvalues that cross the imag¬ 
inary axis are a complex conjugate pair. For smaller 
molecules with the same structure, the first eigenvalue 
to cross the imaginary axis is an isolated real value, 
producing a stationary instability instead of an oscilla¬ 
tory one. It is nevertheless possible that other smaller 
molecules with other structures or with more components 
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FIG. 4. Top: Power Spectrum Density (PSD) of y in ar¬ 
bitrary units for a a = 1.5 and JIa — 0 . 6 , and correspond¬ 
ing to Dc = 5 X 10 “^ (red), Dc = 5 x 10 “^ (green), and 
Dc — 2.5 X 10 “^ (blue). (Inset) Temporal evolution of y, 
which measures the instantaneous deformation of the cluster, 
for JIa = 0.6 and Dc = 5 x 10 “^. Bottom: Power Spectrum 
Density of y in arbitrary units for a a = 1.5 and Dc = 10 “^, 
and corresponding to JIa = 0.7 (red), JIa = 0.6 (green), and 
JiA = 0.5 (blue). 

may present an oscillatory instability. We also note that 
at equilibrium {a a = J^a) the matrix associated with the 
linear dynamics is real and symmetric and therefore all 
eigenvalues are real, and hence oscillations will not be 
possible as expected. 


IV. RUN-AND-TUMBLE MOTION 

An alternative way to exhibit dynamic structure is ob¬ 
tained in molecules that could acquire more than one sta¬ 
ble configuration, dynamically switching between them. 
A good example of this behaviour is demonstrated with 
the AD 3 molecule. In a wide range of parameters, this 
molecule is stabilized with the A colloid located at the 
centre and the three B colloids maintaining contact with 
A, and otherwise free to move in the angular direction as 
sketched in Fig. Depending on the parameters, the B 
colloids can adopt symmetric or asymmetric configura- 
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FIG. 5. Parametrization of the AB^, molecule. The angle 
between the normal vector hi and the x axis is given by (/?, 
whereas 02 and 03 measure the opening angle between con¬ 
secutive B colloids. 


tions resulting in isomers than can self-propel or remain 
Brownian. As in the case of the oscillatory instability, 
for the purpose of analyzing the possible configurations, 
their stability, and the transitions between these config¬ 
urations, we consider the case where the temperature is 
sufficiently low such that we can restrict the dynamics 
by imposing that the B colloids always remain in con¬ 
tact with A (i.e. no fluctuations occur in the relative 
A — B distances). Under this approximation, the equa¬ 
tions of motion for each colloid (without the noise terms) 
are 


Va = {VoUA - Ui)hi + {VoUA - U2)n2 + (Vq/xa - t^3)fi3 


Vbi = —{VottA — Ui)ni + Vo 


ni - n2 
|ni - n 2 |' 


V) 


Vb2 = -(VoSa - U 2 )n 2 + + V^O 


= —{VoaA — Us)ns + Vb 


|n2 - nsl^ 
n3 - ni 
Ins - nip 


■Uo 


ni - ns 
|ni - nsp 
n2 - ni 
|n2 - nip 
ns - n2 
|ns - n2p 


where the normal vectors Ui are defined in Fig. |^and Ui 
are the magnitude of the constraint velocities, which will 
be eliminated using d’Alembert’s principle. 

The motion of the centre of mass of the molecule can 
be readily obtained as 


V — -(Vyi+V^i+VB2+V^s) — (/^A~<^A)(ni+n2+ns), 

(3) 

where it transpires that the molecule propels only un¬ 
der non-equilibrium conditions {a a ^ J^a) and when the 
conformation is not symmetric (lii + n 2 + ris 7 ^ 0). We 
now study under which conditions the configuration is 
not symmetric. 

To derive the motion of the molecule, we apply 
d’Alembert’s principle (see Appendix [A|), considering the 
following generalized coordinates for the relevant degrees 
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FIG. 6. Phase portrait in the 02 — 03 space in the physical 
domain (02 + 03 < Stt). The phoretic charges are JIa = 
1 (top left), /Ia = 0 (top right), and JIa = 1/2 (bottom), 
while a A = 2 in the three cases. Stable hxed points are 
shown in green and saddle points in red. The diagonal solid 
lines indicate configurations with two equal angles. The stable 
conhgurations that correspond to the Y- and T-isomers are 
indicated by arrows. 


of freedom: the x and y coordinates of the central A col¬ 
loid and the three angles 02, and 0s defined in Fig. 

The angle ip is cyclic, due to the rotational invariance, 
implying the conservation law + (02 + ^s + ^(cos pi + 
cos 9^2 + cos (/9s) — y{smpi + sin(/92 + sin(^s) = 0, where 
Pi = p^ P 2 = p ^ (t )2 and (/^s = (/^ + 02 + 0s give the 
orientation of the three B colloids with respect to the x 
axis. The remaining equations of motion can be reduced 
to two coupled equations for 02 and 0s, the bond angles 
of the molecule: 02 = Cc; 2 ( 02 , 0 s) and 0s = (^s( 02 , 0 s)- 
The expressions for UJ 2 and 003 are straightforward to ob¬ 
tain but are quite involved and we do not show them 
explicitly. Figure presents the phase portrait in the 
02 — 0s space where different fixed points appear de¬ 
pending on the values of the phoretic charge JIa, while 
keeping a a fixed (a similar picture is obtained by vary¬ 
ing a A and fixing Jia)^ In Fig. I^top -left, there is only 
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one stable fixed point corresponding to the symmetric Y- 
isomer (02 = 03 = 27r/3), which does not self-propel. In 
Fig. I^top-right, a second stable equilibrium with triple 
degeneration appears, corresponding to the asymmetric 
T-isomer, which self-propels in the direction indicated 
by the double arrow in Fig. [^bottom. Finally, Fig. 
bottom presents the case of bistability between the two 
isomers for an intermediate value of JIa • The bifurcation 
between these phases occurs via the collision of the stable 
fixed points with the saddle points. 

To analyze the transitions between these states in the 
presence of noise we make a further simplification by not¬ 
ing that both the stable fixed points and the saddle points 
lie on the lines where the molecule is partially symmetric, 
with two bond angles being equal. Therefore, we consider 
the reduced dynamics of a single variable 0 = 02 = 03 , 
which is governed by the following equation 

i, , t f~ ~ x(2cos(/) + l)sin(/) 

=UJ{4>) +i = -{otA - MA -- 

3 + cos 20 

sin^ 0 -h sin 0/ (1 — cos 

(3 + cos 20) 

where ^ is the resulting white noise after apply¬ 
ing d’Alembert’s method, with intensity 2I2c/(2 — 
sin0^)cr^. The corresponding Fokker-Planck equation 
(using Stratonovic convention [38]) is 

— - — [ , d ( P \ 

dt 50 ^2 — sin 0^ 50 y ^2 — sin 0^ j 

( 5 ) 

The resulting stationary distributions that is obtained 
by solving equation <§ will depend on uj the Dc- How¬ 
ever, in the low noise limit, the distributions are peaked 
around the stable points, which we can obtain from uj 
only. We find that in the range 0 < JIa < 0.0754 
there is only one stable point at a value of 0 which is 
< 27r/3, and corresponds to the T-isomer. In the range 
0.0754 < JiA < 0.5962 two stable equilibria exist, cor¬ 
responding to the coexistence of the T- and Y-isomers. 
Finally, in the range 0.5962 < JIa the only stable point 
is 0 = 27r/3, namely the Y-isomer. 

In the bistable region, the molecule can exist in the 
two isomeric configurations, making transitions between 
them due to noise. The dynamics is similar to the run- 
and-tumble motion performed by swimming bacteria as 
there is an alternation between run periods performed 
by the T-isomer that are followed by near pauses during 
which the molecule acquires the Y-isomer and undergoes 
rotational diffusion only, leading to the start of a new 
run phase with a random orientation. Here, the new 
orientation is decided by two complementary processes. 
First, in the Y-isomer configuration, there is the rota¬ 
tional diffusion of the cyclic variable cp. Second, when 
the molecule transits to the a new T-isomer, it can adopt 
any of the three degenerate configurations. When the 
rotational diffusion component is small, tumbling occurs 
to three possible discrete new orientations, which is in 


+ 


cos 0/ 




+e, (4) 



x/a 



Vot/a 


FIG. 7. Top: Trajectory of the centre of mass showing the 
run-and-tumble motion of an AB 3 molecule for a a = 2, 
JIa — 0.48, and Dc — 2.5 x 10“^. The open and closed circles 
indicate the beginning and the end of the trajectory, respec¬ 
tively. Bottom: Coarse grained speed of the centre of mass 
for the trajectory in the top panel. The dashed line indicates 
the threshold used to discriminate the run and the tumble 
phases. Five runs and four tumbles are identified as seen in 
the top panel. 


contrast to bacteria that tumble with a continuous dis¬ 
tribution of angles [8]. A sample trajectory of the centre 
of mass is shown in Fig. [^top where the run-and-tumble 
motion is evident, and a video of this motion is presented 
in the supplementary material [39]. The two phases can 
be discriminated by looking at the coarse grained speed, 
which adopts values in the vicinity of 0.6Vb during the 
run phase while it is close to zero in the tumble phase 
(Fig. [^bottom). 

In the bistable region, we have calculated the duration 
of these phases using the Brownian dynamics simulation 
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FIG. 8. Average duration of the run (solid symbols) and 
tumble (open symbols) phases as a function of JIa for fixed 
QA = 2. The noise intensities are Dc = 0.003 (squares) and 
Dc — 0.005 (circles). 


by imposing a threshold value on the coarse grained speed 
(as shown in Fig. [^bottom). The average times are pre¬ 
sented in Fig. The residence times in the run and 
tumble phases exhibit an exponential dependence on the 
value of /iA, while they increase by reducing the noise 
intensity. There is a crossover at JIa ~ 0.51 where both 
phases have similar durations. These behaviours are con¬ 
sistent with what we expect from Kramer’s first passage 
time theory m- 

V. CONCLUDING REMARKS 

We have shown that active colloidal molecules can de¬ 
velop rich dynamical structures, e.g. by having multiple 
stable fixed points or limit cycles in the space of their con¬ 
formations. A key ingredient for obtaining these dynam¬ 
ical states is that the colloidal molecules are not rigidly 
assembled; since they are self-assembled, they are flexible 
and have internal degrees of freedom that can be active 
in the same manner as the global translation or rotation 
modes. 

While we focussed on the simplest examples that ex¬ 
hibit temporal structure, we expect more complex tem¬ 
poral patterns of behaviour to appear for larger and more 
complex molecules. An exemplar movie is presented in 
the supplementary material where random self-assembly 
leads to the formation of a sufficiently large oscillator 
that can break time reversal symmetry and self-propel as 
a whole while beating its oscillatory tail, exhibiting a pat¬ 
tern of motion that is reminiscent of swimming sperma¬ 
tozoa [39] . Further work is needed to study the full range 
of possible structures that can be achieved via nonequi¬ 
librium self-assembly of catalytically active colloids. 

The simple model used for the active colloids serves 
as a proof of concept for self-assembly of molecules with 
dynamic function, namely micro-oscillators and run-and- 
tumble motion. The model has several simplifications. 


notably it neglects hydrodynamic interactions and con¬ 
siders only the far-held phoretic interactions. These and 
other effects such as specihc details of the catalytic re¬ 
action kinetics should be included in order to have a re¬ 
alistic description of the system that allows us to make 
quantitative predictions. Therefore, in practice, it is not 
guaranteed whether the particular examples presented in 
this manuscript will perform their time dependent func¬ 
tion when all these effects are included. However, these 
details will only make quantitative changes, which means 
that similar conhgurations will exhibit the reported func¬ 
tionality, because those effects do not alter the main un¬ 
derlying mechanism for dynamic function, which is the 
action-reaction symmetry breaking that is present in the 
phoretic interactions. 


Appendix A: D’Alembert’s principle in overdamped 
dynamics 

Consider an ensemble of particles with overdamped dy¬ 
namics 

ri=Ff; i = (Al) 

where is the total force acting on each par¬ 

ticle and F^ are constraint forces, where the mobilities 
have been set to one for simplicity. In the case of col¬ 
loidal particles, these constraint forces keep the relative 
distances of the particles in contact fixed. 

Generalized coordinates g/., k = 1,2,...,n < A, are 
defined satisfying the constraints such that the particle 
positions can be uniquely calculated in terms of them; i.e. 
Ti = r^(gi,..., g^). With these generalized coordinates 
virtual displacements Svi can be constructed such that 
they satisfy, by construction, • FJ • dvi = 0. Therefore, 

(Fj - ii) ■ 6ri = 0. (A2) 

i 

Using dvi = the previous expression can be 

written as 

^ (^Qk - Sqk = 0, (A3) 

where 

i 

are the generalized forces, which do not depend on the 
constraint forces, and 

T = Y1 (A5) 

i 

plays a role analogous to the kinetic energy. 





As the virtual displacements in the generalized coordi¬ 
nates are independent, the resulting equations of motion 
in the generalized coordinates are 

^T = Qu] A:= l,2,...,n. (A6) 

dqk 
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